Repulsively-bound exciton-biexciton states in high-spin fermions in optical lattices 
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We show that the interplay between spin-changing collisions and quadratic Zeeman coupling 
provides a novel mechanism for the formation of repulsively bound composites in high-spin fermions, 
which we illustrate by considering spin flips in an initially polarized hard-core ID Mott insulator 
of spin-3/2 fermions. We show that after the flips the dynamics is characterized by the creation 
of two types of exciton-biexciton composites. We analyze the conditions for the existence of these 
bound states, and discuss their intriguing properties. In particular we show that the effective mass 
and stability of the composites depends non-trivially on spin-changing collisions, on the quadratic 
Zeeman effect and on the initial exciton localization. Finally, we show that the composites may 
remain stable against inelastic collisions, opening the possibility of novel quantum composite phases. 

PACS numbers: 



I. INTRODUCTION 

Ultra cold gases in optical lattices offer an extraordi- 
nary controllable environment for the analysis of many- 
body phenomena, as highlighted by the realization of 
the Mott insulator to superfluid phase transition in Bose 
gases Qj . and more recently by the observation of the 
Mott insulator to metal transition in two component 
Fermi systems HM- I n addition to a wealth of possible 
quantum phases and due to its discrete nature, op- 

tical lattices offer as well novel possibilities for composite 
formation, including local on-site dimer (or n-mer) for- 
mation in attractive gases @ and fermionic composites in 
Bose-Fermi mixtures [7] . An even more striking effect of 
the lattice discreteness is given by the recent observation 
of repulsively bound atomic pairs (doublons) at a given 
lattice site [g, which are dynamically stable in the 
absence of dissipation. Note that contrary to standard 
condensed-matter systems, where e.g. interactions with 
phonons would lead to a rapid dissipation of the dimers, 
optical lattices provide to a large extent a dissipation- 
free environment. Indeed experiments on optical lattices 
provide an extraordinary scenario for the study of the 
intriguing physics of metastable bound states and other 
far-of-equilibrium phenomena 

Spinor gases, formed by atoms with various available 
Zeeman substates, present a very rich physics due to 
the interplay between internal and external degrees of 
freedom Interestingly, interatomic interactions lead 
to spin-changing processes in which population is trans- 
ferred between Zeeman sublevels. The corresponding 
spinor dynamics has attracted a large interest, m ostly 
in the realm of spinor Bose-Einstein condensates [Il| - [l3j . 
Spinor lattice gases offer fascinating novel physics, most 
relevantly in what concerns quantum magnetism, includ- 
ing antiferromagnetic order in spin- 1/2 fermions l4 and 
even more intriguing phases for higher spins (TBI - 20| . 

Repulsively on-site bound pairs in spin- 1/2 fermions 
are formed since the large interaction energy of the dou- 
blon cannot be accommodated by the maximal kinetic 
energy for two atoms in the lowest band (proportional to 
the hopping energy) @, [2l| - |23| . In this paper we show 
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FIG. 1: (Color online) (a) Exciton of a 3/2 particle and a 
—3/2 hole; (b) Biexciton of by two m = ±1/2 particles and 
two m — —3/2 holes; (c) Broken biexciton due to C2 processes 
(see text). 

that a completely different mechanism, based on the in- 
terplay between spin-changing and quadratic Zeeman ef- 
fect (QZE), may sustained novel types of composites for 
repulsive high-spin lattice fermions. We illustrate this 
physics for the case of hard core ID spin-3/2 fermions 
in a Mott insulator initially polarized in m = —3/2. We 
show that spin- flips into m — 3/2 may lead to two types 
of composites formed by an exciton (m = 3/2 particle 
and m = —3/2 hole) and an antisymmetric biexciton 
(m = ±1/2 particles and two m = —3/2 holes) (Fig. [IJ. 
We show that the dynamics and stability of the exciton- 
biexciton composites exhibits a non-trivial dependence 
with spin-changing, QZE and center-of-mass momentum. 
Finally, we show that inelastic composite-composite in- 
teractions may be very inefficient, opening exciting pos- 
sibilities for many-body composite gases. 

The structure of the paper is as follows. In Sec. [TTI we 
introduce the spin-3/2 system under consideration as well 
as the possible exciton, and biexciton excitations occur- 
ring after an individual spin flip. Sec. Mil introduces the 
idea of repulsively bound exciton-biexciton pairs, with a 
particular emphasis in their band-like dispersion and the 
conditions for their existence. Sec. IIVI is devoted to the 
intriguing dynamics of the repulsively bound pairs, and 
its dependence on the spatial derealization of the ini- 
tial spin-flip excitations. In Sec. [V] we analyze the case 
of multiple spin-flip excitations. Finally, in Sec. I VII we 
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summarize our conclusions. 



II. MODEL 



A. Hamiltonian 



stant c 3 /2 (ci/2) characterizes the super-exchange be- 
tween neighboring sites with ±3/2 (±1/2), whereas c 2 
does the same for sites with to and m' such that \m\ 7^ 
\m'\. Finally, c sc characterizes the super-exchange lead- 
ing to spin-changing from ±1/2 to ±3/2 or viceversa. 



We consider spin-3/2 fermions (to = ±3/2, ±1/2) in 
a deep ID lattice, such that at low filling only the low- 
est band is relevant. In this regime, the physics is given 
by the interplay between inter-site hopping (character- 
ized by the constant t), QZE (given by a constant q, 
which may be externally controllable), and s-wave col- 
lisions. The latter may occur in two different channels 
with total spin F = and 2, characterized by the cou- 
pling constants gF = Aith^ap /M, with of the scattering 
length and M the atomic mass. Although typically ao 
and CL2 are similar, they may be varied by microwave 
dressing [24j or optical Feshbach resonances [2a ], Below 
we use G = (g a +g 2 )/2, and g = (g 2 - 9o)/{92 +9o)- The 
s-wave interactions preserve magnetization (and hence 
the linear Zeeman effect is irrelevant), but may induce 
spin-changing collisions (characterized by g), which re- 
distribute the populations at different components. For 
G>i we may consider maximally one fermion per site 
(hard-core case). For a sufficiently large chemical poten- 
tial the system enters into the Mott insulator with one 
fermion per site. The spin physics within the Mott insu- 
lator is given by an effective Hamiltonian: 



H = qm 2 h m .i + 2J H UJ , 
where (■ • • ) denotes nearest neighbors and 
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B. Single-flip excitations 



We consider a system initially polarized into to 



-3/2, \1> BG ) =Uii> 



-3/2, 



\vac) 



where \vac) is the vac- 



uum. \iPbg) is stable since collisions preserve magne- 
tization. We are interested in excitations of the form 
\m,j) = ^ ? t_3/2j'0-3/2,j#BG), created by spin flips. 
Note that these excitations have actually an exciton- 
like character, since they are formed by a particle (in 
to ^ -3/2) and a hole (in m = -3/2) (Fig. [Ik). 

A single spin-flip into to = ±1/2 is not accompanied 
by spin-changing, and hence lead to isolated ±1/2 ex- 
citons which hop with an effective hopping C2. Spin- 
changing collisions lead to a much less trivial dynamics 
when spins are flipped into m = 3/2. After the spin-flip 
the 3/2 exciton may hop to the nearest neighbor via c 3 / 2 
super-exchange (second line in Eq. (0)). In absence of 
spin-changing {g — 0) the 3/2 exciton propagates freely 
with an energy E 3 / 2 (k) = 2c 3 / 2 (l — cos k), where k is the 
center-of-mass momentum of the 3/2 exciton. 

Spin-changing modifies this simple picture, since the 
to = 3/2 spin may interact via c sc super-exchange with a 
neighboring —3/2 spin leading to two neighboring spins 
with to = ±1/2, i.e. a biexciton of the form |m, j; —m,j+ 

1) =) = V4, i ^ rn , i+1 V>-3/2, J - + lV'-3/2, J #BG> (Fig. lb). 

After spin-changing, one of the newly created excitons, 
say to = 1/2, may move apart from the neighboring 
to = —1/2 one by a c 2 superexchange (Fig. [TJ;). Fur- 
ther C2 processes bring the pair even further away from 
each other. Since an isolated ±1/2 spin surrounded by 
—3/2 neighbors remains stable against spin-changing, 
one naively expects all 3/2 excitons to eventually dis- 
solve into independent ±1/2 ones, which would freely 
move away with a hopping c 2 . As shown below, contrary 
to the naive expectation, the 3/2 excitons may become 
very robust. 



describes the effects of the interatomic interactions, with 
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c sc = t 2 sin 2(j> ^2 rj , 

-t 2 /{g + G), 
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(3) 
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where rj ± = ±1, A± = G - hq/2 ± [4q 2 + g 2 G 2 ] 1/2 , 
and tan0 = ([4q 2 + g 2 G 2 ] 1 / 2 + 2q) / gG [26]. The con- 



III. EXCITON-BIEXCITON REPULSIVELY 
BOUND STATES 

A. Band structure 

We introduce the Fourier transformed states [27j 
\4>f{k))^^Tti^ Hl+3/2) \ T 1/2, i; ±1/2,/ + j) and 
\ip3/2{k)) = -^j- z\jf=i e ikl \3/2, j), where L is the number 
of sites and k the center of mass quasi- momentum either 
of the 3/2 exciton or of the ±1/2 biexciton. Note that 
the biexciton states \4>, (k)} contain not only biexcitons 
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FIG. 2: Band spectrum for single-flip excitations as a function 
the quasimomentum k, for q = —4 and g — 0.8. The circles 
indicate the results of Eq. QQ) . 



of neighboring ±1/2 (i.e. j = 1, as in Fig. [T]d), but also 
broken biexcitons (i.e. j > 1, as in Fig.Ht). This is true 
even for the bound solutions discussed before, which may 
contain as well a contribution of broken biexciton pairs. 

It is particularly useful to introduce the symmet- 
ric and antisymmetric states \4>j' A (k)) — (\<t>t(k)) ^ 
\4>j{k)))/y/2. Note that the Hamiltonian © do not 
couple symmetric and antisymmetric states, and hence 
the Hamiltonian in fc-space splits into H = J2 k (Hs{k) + 
H A {k)). The physics of the symmetric biexcitons is given 
by 

H s (k) = E ls \cf>!(k))(c(>f(k)\+E B J2\^ik)){^(k)\ 

- n 2 (k)J2(\^(k))(^ +1 (k)\+H.c), (7) 

i>i 

with Eis = 2c 2 - 4q, E B = 4c 2 - 4g and fi 2 (fc) = 
2c2 cos(fc/2), whereas the coupling between 3/2 excitons 
and antisymmetric ±1/2 biexcitons is given by 

H A (k) = E 1A \<t>?(k))(<fi?(k)\+E B J2\<t>f(k))(<t>f(k)\ 

+ E 3/2 (k)\ij 3/2 (k))(<ip 3/2 (k)\ 

+ O sc (fe)(# 3 / 2 (fc))^l(fc)|+-ff.C.) 

- ^ 2 (fc)^(|^(fc))(^ +1 (fc)| +J ff. C .), (8) 

where = 2c 2 + 2c 1/2 - 4g, E 3/2 (k) = 2c 3/2 (l - cos fc), 
and fi sc (£;) = v^Csc sin(fc/2). 

Figure [5] shows a typical spectrum (which resem- 
bles that of self-bound repulsive pairs of spin-1/2 
fermions Q). Unbound ±1/2 biexcitons form a con- 
tinuous band, due to their free relative motion with 
an effective hopping f2 2 (fc), which leads to a dispersion 
E(k,k r ) = Eb + 2tt 2 (k) cos k r , where k r is the relative 
momentum of the pair. We observe also the appear- 
ance of up to three bound states (resembling the case 
of repulsiv e sp in-1/2 fermions with nearest neighbor in- 
teractions 27]). One of the bound states (upper state 



in Fig. [2]), with energy E(k) = E B - 2c 2 (l + cos 2 (fc/2)), 
belongs to the symmetric manifold, being hence decou- 
pled from the 3/2 exciton, and thus we do not consider 
it further. 



B. Pairing of excitons and antisymmetric 
biexcitons 

The other two bound states result from the coupling 
of antisymmetric ±1/2 biexcitons and 3/2 excitons. In 
order to understand these states, we introduce the ansatz 

\m) = cos^^^ir^-^'- 1 )/^^^)) 



i sin ¥#3/2 (fe)), 



(9) 



where 7 is the localization length. Imposing E\i[i(k)) = 
H A (k)\ijj(k)) results in a transcendent equation for 7: 



E B + 2ft 2 cosh(l/ 7 ) = 
~ \E 1A + Q 2 e 1 ^ + E 3/2 



±^[E 1A + n 2 eMt ~E 3/2 y 



(10) 



where the relation between energy and 7 is given by 

E{k) - E B = 2fi 2 (fc)cosh(l/ 7 (fc)). (11) 



In Fig. [2] we indicate the results obtained from Eqs. (fT0|) 
and (fTTj). which are in excellent agreement to those di- 
rectly obtained from the diagonalization of ©. Eq. (fTTj) 
shows that the separation between the bound state and 
the lowest boundary E B + 2fl 2 (k) of the sea of unbound 
±1/2 pairs is crucial for its localization, and hence while 
the lower bound state is tightly bound, the upper one is 
much more loose. In particular, the upper bound state 
unbinds (7 diverges) if its energy approaches E B + 257 2 . 
Hence the existence of the bound exciton-biexciton states 
depends non-trivially on QZE, spin-changing collisions, 
and exciton momentum. Fig. [3] shows those values of g 
and q for which one or two bound states are expected for 
k = 0. 



IV. DYNAMICAL PROPERTIES OF THE 
EXCITON-BIEXCITON BOUND PAIRS 



The mixing angle (p fulfills 



tan ip(k) = 



1 



^SC (k) 



E(k) - E 1A - n 2 (k)e 



-1/7W 



(12) 

The contribution of the 3/2 exciton to the bound states 
(sin 2 ip) is particularly important to understand the dy- 
namics following a spin- flip into m = 3/2. In particular, 
at k — (when the center of mass of the spin-flip exci- 
tation is maximally delocalized), fl sc — and hence the 
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a linear superposition of both bound states (and possi- 
bly a slight contribution of broken biexcitons). Since the 
bound states have different group velocities of opposite 
sign, the dynamics after the spin-flip is hence charac- 
terized by the appearance of two wavepackets moving 
in opposite directions formed by the two different types 
of exciton-biexciton composites (Fig. [SJ middle, which 
should be compared with the single wavepacket spread- 
ing for k — in Fig. O top). 

The situation is markedly different if the initial spin- 
flip is strongly localized, where all k in the Brillouin zone 
contribute. As a consequence, the initial 3/2 excitation 
results in an admixture of both bound states, which re- 
sult, as shown in Fig. [5] (bottom), to a completely differ- 
ent wavepacket dynamics as that for initially extended 
spin- flips (Fig. El top). Note that, in addition, a localized 
initial spin-flip may become considerably more unstable 
against unbinding, especially if the upper antisymmetric- 
3/2 bound state unbinds for a given range of k values. In 
that case a significant part of the initial exciton excita- 
tion would dissolve into unbound biexciton pairs. 



V. MULTIPLE SPIN-FLIP EXCITATIONS 



FIG. 4: Time evolution of the width (j 2 ) of the 3/2 exciton 
wavepacket for t — 0.1G, g = 0.8 and qt 2 /G = —4 with (dot- 
ted) and without (solid) spin-changing collisions, for an initial 
Gaussian wavepacket with a width of 10 sites. 



upper bound state is a pure antisymmetric ±1/2 biexci- 
ton with E = E 1A + 4c%/E 1A and 7 = 1/ ln(E 1A /2c 2 ), 
whereas the lower bound state is a pure m = 3/2 exciton 
with energy £^3/2(0). On the contrary, when the center of 
mass of the spin-flip excitation is strongly localized there 
is a strong coupling between the exciton and the anti- 
symmetric biexciton which leads to a strongly distorted 
dynamics, as shown below. 

If the initial spin-flip is delocalized in a spatial region 
much larger than the inter-site spacing the momentum 
distribution is strongly peaked at k = 0. Hence after 
the spin-flip there is no significant transfer into ±1/2 
biexcitons, since the lower bound state is basically a 3/2 
exciton. Without spin-changing the exciton tunnels with 
hopping C3/2, and hence at k = the effective exciton 
mass is m*o = l/2 c 3/2- In spite of the absence of any 
significant biexciton admixture, spin-changing strongly 
modifies the effective exciton mass into: 



Note that m* differs significantly from 771*0 (it may 
even invert its sign), radically modifying the exciton 
wavepacket dynamics, as illustrated in Fig. [4] 

Note that if the created initial exciton has a sharply 
defined but finite center-of-mass momentum k , then the 
3/2 exciton is not any more a single bound state but 



Up to this point we have considered the dynamics af- 
ter single spin-flips. Although single spin flips may be 
created using state of the art techniques for single-site 
addressing [28|, [29| . in typical experiments more than 
one spin-flip will be induced. As a result more than one 
composite will be simultaneously created, hence opening 
the possibility of inelastic composite-composite collisions 
which may induce a decay channel for the bound com- 
posites into broken biexciton pairs. The various initial 
spin-flip excitations result in a complicated many-body 
non-equilibrium dynamics for the created bound compos- 
ites, which we have studied by means of time-dependent 
Matrix-Product-State calculations [3(| for up to 5 spin- 
flips delocalized in an initial extension of 10 sites (for a 
total of 30 sites). 

Figure [B] shows our results for the case of qG/t 2 = —4 
and g = 0.6. For this particular case the bound states are 
well below the sea of broken pairs. Hence, as discussed 
above, it is expected that single spin-flip excitations do 
not lead to a significant population of unbound solutions. 
Two main features can be observed in Fig. [51 On one side, 
there are marked Rabi like oscillations of the 3/2-exciton 
population due to the presence of the two bound states. 
The oscillation amplitude is small since, as discussed in 
the previous section, the lowest bound state is to a large 
extent the 3/2-exciton due to the initial delocalization of 
the spin-flip. Note that the Rabi oscillations are how- 
ever damped due to dephasing, since different momen- 
tum components close to fco = have a slightly different 
Rabi frequency. Note that a more localized initial spin- 
flip would result in a stronger frequency beating, due to 
the participation of momenta over all the Brillouin zone. 
Note also that the presence of additional spin-flips leads 
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FIG. 5: Time evolution of the magnetization rm = 
12m m (^™.,i) for an initial Gaussian wavepacket for the flip ex- 
citation with width of 10 sites and central momentum ko = 
(top), and fco = n/2 (middle), and for a localized initial flip 
with width of 1 site and fco = (bottom) 

to an even stronger damping of the Rabi oscillations. 

On the other side, note that composite-composite col- 
lisions basically leave unaffected the total population in 
the bound states, which for this particular case basi- 
cally reduces to the population of 3/2 excitons (Fig. [1^) 
and unbroken ±1/2 biexcitons (Fig. [TJd). Hence, deeply 



FIG. 6: (Color online) Time evolution of the bound popula- 
tion (top) and of the population of 3/2 excitons (bottom) for 
t = 0.1G, g = 0.6 and qt 2 /G = -4, for the case of 1 (dotted), 
3 (dashed) and 5 (solid) spin flips delocalized in an initial 
square wavepacket 10 sites wide. 

bound states may be created, which do not undergo sig- 
nificant inelastic unbinding due to inelastic collisions. 
Thus, multiple spin-flips lead to an out-of-equilibrium 
(but highly metastable) many-body state, which may un- 
derstood as a basically stable gas of bound composites, 
opening interesting perspectives for quantum composite 
gases. For the particular case considered in Fig. [6] we 
expect a quantum gas of almost pure 3/2-excitons with 
"biexciton dressed" dynamics. Note that although the 
3/2-excitons present an effective attractive nearest neigh- 
bor interaction (— 2c 3 / 2 ), this interaction becomes irrel- 
evant at low momenta k — ¥ compared to the infinite 
on-site repulsion due to Pauli exclusion. Hence the low k 
properties of ID 3/2 excitons will be as those of a Tonks- 
Girardeau gas. 



VI. CONCLUSIONS 

In summary, spin flips in a polarized repulsive high- 
spin Fermi gas may lead to the formation of novel types 
of bound composites, which we illustrated for the case 
of ID spin-3/2 hard-core fermions in the Mott phase. In 
that case the composites are formed by a 3/2 exciton- 
like excitation and an antisymmetric ±1/2 biexciton-likc 
one. Intriguing dynamics and stability properties of the 
composites result from a non trivial interplay between 
spin-changing, QZE and exciton momentum. The stabil- 
ity of the exciton gas against inelastic interactions opens 
exciting possibilities for the creation of intricate novel 
quantum composite phases. 



Acknowledgments 

We thank the DFG for support (Center of Excellence 
QUEST). 



6 



[1] M. Greiner et ai, Nature (London) 415, 39 (2002). 

[2] R. Jordens et ai, Nature (London) 455, 204 (2008). 

[3] U. Schneider et al, Science 322, 1520 (2008). 

[4] M. Lewenstein et al, Adv. Phys. 56, 243 (2006). 

[5] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[6] A. Rapp, G.Zarand, C. Honerkamp, and W. Hofstetter, 

Phys. Rev. Lett. 98, 160405 (2007). 
[7] M. Lewenstein et al, Phys. Rev. Lett. 92, 050401 (2004). 
[8] K. Winkler et al.,, Nature (London) 441, 853 (2006). 
[9] N. Strohmaier et al, Phys. Rev. Lett. 104, 080401 (2010) 
[10] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 
[11] H. Schmaljohann et al., Phys. Rev. Lett. 92, 040402 
(2004). 

[12] L. E. Sadler et al, Nature 443, 312 (2006). 
[13] C. Klempt et al, Phys. Rev. Lett. 104, 195303 (2010). 
[14] F. Werner et al, Phys. Rev. Lett. 95, 056401 (2005). 
[15] C. Honerkamp and W. Hofstetter, Phys. Rev. Lett. 92, 

170403 (2004). 
[16] C. Wu, Phys. Rev. Lett. 95, 266404 (2005). 
[17] P. Lecheminant et al, Phys. Rev. Lett 95, 240402 (2005); 

S. Capponi et al, Phys. Rev. B 75, 100503(R) (2007); G. 



Roux et at, Eur. Phys. J. B 68, 293 (2009). 
[18] Gorshkov et al, Nature Physics 6, 289, (2010). 
[19] M. A. Cazalilla et al, New J. Phys. 11, 103033, (2009). 
[20] R. Jordens et al, Phys. Rev. Lett. 104, 180401 (2010). 
[21] D. Petrosyan et al, Phys. Rev. A 76, 033606 (2007) 
[22] F. Alzetto and X. Leyronas, Phys. Rev. A 81, 043604 

(2010) 

[23] Y.-M. Wang and J.-Q. Liang, Phys. Rev. A 81, 045601 
(2010) 

[24] D. J. Papoular, G. V. Shlyapnikov, and J. Dalibard, 
Phys. Rev. A 81, 041603(R) (2010). 

[25] P. O. Fedichev et al, Phys. Rev. Lett. 77, 2913 (1996). 

[26] A divergence occurs when q = ±G(g 2 — l)/4, where per- 
turbation theory breaks down, since the Mott insulator 
energy equals that with a particle-hole excitation. We 
avoid this divergent case. 

[27] J.-P. Nguenang and S. Flach et al, Phys. Rev. A 80, 
015601 (2009). 

[28] W. S. Bakr et al, Nature 462, 74 (2009). 

[29] J. F. Sherson et ai, Nature 467, 68 (2010). 

[30] F. Verstraete et al, PRL 93, 227205 (2004). 



